Fermionic trimers in spin-dependent optical lattices 



O 

(N 



O 



X 



Giuliano Orso 1 , Evgeni Burovski 2 , and Thierry Jolicoeur 2 
1 Laboratoire Materiaux et Phenomenes Quantiques, 
Universite Paris Diderot-Paris 7 and CNRS, UMR 7162, 75205 Pans Cedex 13, France and 
2 LPTMS, CNRS and Universite Paris-Sud, UMR8626, Bat. 100, 91405 Orsay, France 

(Dated: October 15th, 2010) 

We investigate the formation of three-body bound states (trimers) in two-component Fermi gases 
confined in one dimensional optical lattice with spin-dependent tunneling rates. The binding energy 
and the effective mass of the trimer are obtained from the solution of the Mattis integral equation 
generalized to the case of unequal Bloch masses. We show that this equation admits multiple 
solutions corresponding to excited bound states, which are only stable for large mass asymmetry. 



I. INTRODUCTION 



The advent of Feshbach resonances and optical lattices has caused a major revolution in the field of ultra-cold atoms 1 . 
Few-body physics is currently playing an important and intriguing role. First, bound states of few interacting atoms 
c/2 . can be produced and studied experimentally in a controlled way, providing a direct test to fundamental quantum- 
mechanics. Recent examples are the observation of repulsively bound pairs of bosonic atoms in an optical lattice 2 - 
and the detection of an Efimov-like trimer 3 of 133 Cs atoms. 

Second, optical lattices can considerably affect the many-body scenario by modifying the two-body scattering 
properties of interacting atoms^. The connection between few- and many-body physics is particularly interesting in 
one dimensional (ID) systems, where the strongly interacting regime typically occurs at low density. After the first 
q-i observation^ of confinement-induced molecules in attractive Fermi gases, experimentalists are now addressing the rich 
! t . many-body scenario predicted in these systems. 

In this context, a very recent experiment performed at Rice University investigated the properties of a one dimen- 
sional trapped Fermi gases with attractive interaction and a finite spin polarization, verifying the two-shell structure 
of the density profiles predicted in Refs.— ^ on the basis of Bethe-Ansatz calculations. In particular the partially 
polarized core of the gas is expected to be a supcrfluid of the Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) typ o 10 ' 11 , as 
confirmed by exact numerical simulations^. 

Motivated by the strong interest in this field, in Refsi 13 ' 14 we investigated the ground state properties of a two- 
component Fermi gas confined in a ID lattice with spin-dependent tunneling rates. The system is described by the 
asymmetric Fermi-Hubbard model — ~— : 



H = - J~] t a (c\ a c i+1 ^ + h.c)j + U ^ i (!) 



(N 
> 
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. where U < is the on-site attraction and t a are the spin-dependent tunneling rates. Here Ci„ annihilates a fermion 
with spin a at site i and hi a is the local density. The Hamiltonian (JXJ) can also describe mixtures of two different 
atomic species, like mixtures of 6 Li and 40 K near a heteronuclear Feshbach resonance^—. 

For equal tunneling rates, t± — tf, the exact Bethe Ansatz solution of the model ([1} shows that n-body bound 
I states with n > 2 are generally forbidden 2 ^. In Ref. 1 - we showed that in the presence of unequal tunneling rates, 
• • ■ there is formation of trimers made of two heavy (4.) fermions and one light (f) fermion. These states are reminescent 
of the trions recently observed 2 ^ in semiconductors. Then a DMRG calculation shows that trimers are responsible for 
the appearance of a new gapped phase for finite and commensurate densities (n± — 2n-|-), which is characterized by 
$H ' exponential suppression of both single-particle and superfluid FFLO correlations. 

The purpose of this paper is to present a more thorough derivation of the three-body calculations outlined in RefJ^. 
We also discuss new and unexpected results for the excited bound states. The three-body problem of interacting atoms 
with unequal mass has already attracted a lot of attention 2 ^—. We will make contact with the corresponding results 
for the continuum model obtained in Ref^. 

The article is organized as follows. In section [TT] we present a self-contained derivation of the Mattis integral 
equation 3 -! for the three-body problem, generalized to the case of unequal masses. In section IIIII we calculate the 
ground state solution of such equation, corresponding to the trimer state with the highest binding energy. In section 
II VI we investigate the excited bound states solutions of the integral equation. Finally we give our conclusions in section 

E) 

Before continuing, we would like to mention that three-body bound states (though of different nature) also occur 
in the ID Bose-Hubbard model, as recently investigated in Re L 32 ? 33 , as well as in Fermi gases with three or more spin 
components^—. 
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II. INTEGRAL EQUATION FOR TRIMERS 

In this section we consider two .J,- fermions interacting with one f-fermion, as described by the Hamiltonian ([l} . We 
map the Schrodinger equation of the three particles into an integral equation which we then solve both analytically 
and numerically. The Schrodinger equation in momentum space takes the form 

(e^k 1 ) + e i {k 2 ) + e^(k 3 )-E)^{k l ,k2,k 3 ) + U J ^{p,k 2 ,kx + k 3 - p) + U J ^{k uP ,k 2 + k 3 - p) = 0, (2) 

where e a (k) — 1t a {\ — cosfc) the energy dispersions of the two components and the integration over quasi-momenta 
is restricted to [— tt, it]. We see from Eq.© that the total quasi- momentum P = k\ + k 2 + k 3 is a conserved quantity 
associated to the discrete traslational invariance of the lattice model^. Introducing the function 



tf>(p,k,P-p-k), (3) 



and taking into account that the wavefunction is antisymmetric under exchange of the two J. fermions, ip(ki,p,P — 
ki — p) = —ip{p, k\,P — k\ —p), we can rewrite Eq.@ as 

, , x A(k 2 ,P) -A(k!,P) 

ei(ki) + ei{k 2 ) + e t (fc 3 ) - E 

which can be seen as a self consistent equation for the function A(k, P). Inserting Eq.@ into Eq.Q we obtain 

Mk,P) = -U T f A ^p-^f , (5) 
J -ix 2?r £(k,q,P)-E 

where E{k,q 1 P) = t^(k) + ei(q) + e^(P — k — q) is the total energy dispersion. Bringing the term proportional to 
A(k,P) in Eq.© to the right hand side, we get 

A(A, P)(l + UUk, P)) = U £ Z ei^-E ' < 6 » 
where the integral I E (k,P) is defined by 

I B {k,P)= I ^ m * = ■ 
J 2n t{k,p,P) - E 

By setting z = e lp and using the residue theorem for complex functions, we find 



(7) 



I E (k, P) = , 1 = ■ (8) 

fa(k) + 2t x + 2 - Ef - 4 - 4tJ - 8*4 cos(P - k) 

The function I E (k,P) appears already in the solution of the two-body problem in a lattice 3 ^. In particular the 
condition 1 + UlE(k ~ 0, P = 0) — yields the binding energy i?p air = —E of a pair of up and down fermions. From 
Eq.© one finds E b pail = -2(1 + t ; ) + ^JU 2 + 4(1 + t^. 

Introducing R E (k,P) = (1 + UI E (k,P)) 1/2 and the function f(k,P) = A(k, P)R E (k, P), Eq.© takes the form of 
a homogeneous integral equation — 



Uf(q,P) 

R E {k,P)R E (q,P)[£(k,q,P)-E] ' 



(9) 



whose solution yields the energy E of the three-body system. We are interested here on bound states solutions 
corresponding to energy E < —E paiT . 

It is important to notice that for zero total quasi- momentum (P — 0), the function f(k,P = 0) in Eq.Q must 
be odd, namely /(— k, 0) = — f(k, 0). This comes from the fact that in this limit the function ([3]) coincide with the 
Fourier transform of the real space wave-function when two interacting particles are at the same lattice site: 

lKn x ,n a , m) = f ^- e ^"^A(k, P = 0), (10) 
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FIG. 1: Binding energy of the ground state trimer as a function of the tunneling rate tj, for different values of the interaction 
U = — 2(bottom curve), —4, —8. The asymptotic curve [see Eq, (|18|l ] in the strong coupling limit U = — oo is shown with the 
dashed line. 



where ip is a zero quasi-momentum state. Since ip(ni, n%, tlx) = —ip{ri2, n\, rii), we see from Eq. (jlOp that A{k, P = 0) 
is an odd function of the quasi-momentum. Since Re{— k, P = 0) = Rs(k,P = 0), we conclude that the function 
f(k, 0) is also odd. This property will be used systematically below to obtain our analytical results. 

In the following we first discuss the properties of the three-body bound state with the lowest energy E (or, equiv- 
alently, with the largest binding energy). Excited bound states solutions will be discussed next. To simplify the 
notation, from now on, we fix the energy scale by set if = 1 in Eq.©. 



III. GROUND STATE SOLUTION 



A. Binding energy 



Equation ^ can be considered as an eigenvalue problem Ke • f = Af , where the energy E is fixed by the constraint 
A = 1. We solve this equation numerically for zero quasi-momentum P = 0. The binding energy E^ T of the trimer is 
related to the total energy E by — E = E^ air + E^.. In Fig. [I] we plot the binding energy of the trimer as a function 
of the mass asymmetry t± for increasing values of the attraction U . 

We see that E% T vanishes at the symmetric point fj, = 1 for any values of the interaction U, in agreement with 
the Bethe Ansatz solution^. For t± > 1 no bound state solution has been found. This result can be understood by 
noticing that when the two heavy particles approach each other, the light fcrmion can hop between the two without 
loss of potential energy. Therefore, if ij, < 1, the energy gain to delocalize the light particle overcompensates the 
energy cost to localize the heavy fermions, and the trimer state is bound. 

As the mass asymmetry increases the binding energy also increases until it saturates at t± — > 0, where the effective 
mass of the heavy fermions becomes infinite. In this limit the function ((8|) reduces to a constant -Tb(^) = 1/ \J E{E — 4), 
implying that Re{q) = (1 + Z7 j \J E(E — 4)) 1 / 2 = Re in Eq. (|9]) is also constant. By changing the integration variable 
to q' = q + k and omitting the prime index, the latter takes the form 



J[L] = w 



E 



dq f(q - k) 



_ n 2tt 2(1 - cosg) - E 



(11) 



Taking into account that f(k) is an odd function, we can write the general solution as f(k) = Xm=i ° n sin(nq'). 
Substituting this into Eq. pTj) , we find that the different harmonics decouple, implying that the solution is of the form 
fn{k) = sin(nfc). By inserting this into Eq. (fTTj) , the latter reduces to 1 = y n (E), where 



U 

K 



E 



COs(nq) 



K dq 

„ 2^ 2(1 - cosg) - E 



(12) 
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FIG. 2: (Color online) Energy dependence of the function y n (E) [see Eq. (112[l ] plotted for U = —2 and for different values of the 
index n = l(top curve), 2, 3, 4, 5. The dashed vertical line corresponds to the threshold energy — £^ alr — 2 — \/U 2 + 4 = —0.828 
for the pair bound state for ti = 0. The energy of the bound state is found from the solution of the equation 1 = y n (E). The 
ground state corresponds to n = 1, where E = —U 2 /(l — U) = —1.333 [see Eq. (|14|l ]. 



The function y n {E) is plotted in Figf2]for U = —2 and for increasing values of n, starting from n — 1 (top curve). 
The dashed vertical line corresponds to E = — £p ail . = 2 — \fU 2 + 4. We see that there is an infinite number of three- 
body bound state solutions that for large n accumulates near the pair energy. The ground state solution corresponds 
to n — 1, namely f(k) = sinfc. By using the formula J* cosq/(a — cos q)dq = 2Tx{aj\Ja 2 — 1 — 1) valid for a > 1, as 
well as the explicit expression for R E , from Eq.(fT2|) we obtain 



j _ -U (2-E-y/E(E-4) 
2 y/E(E - 4 + U ' 



yielding E = —U /(l — U). Therefore the trimer binding energy at t± = is given by- 

U 2 



ElM = 0) = +2 - VU* +4 , (14) 

in agreement with our numerical results in FigJTJ In this limit the problem is very simple since we are dealing with 
two fixed scatterer and only one mobile particle : this point of view will be used in detail in section 4. 

Let us now discuss the dependence of the binding energy on the interaction U . Clearly, the binding energy increases 
as U increases. However, differently from the pair binding energy, which diverges for infinite attraction, E^ T saturates 
to a finite value shown in Fig. [T]with dashed line. In this strong coupling regime, corresponding to \E\ ~ \U\ 3> 1, we 
can use the expansion (£ — E)^ 1 ~ —l/E — E/E 2 in Eq.(j9]). The first term gives no contribution due to symmetry 
considerations, whereas the second term yields 

r dq Ucos(q + k)f(q) 
m -J_^ E 2 R E (k)R E (q) • (15) 

By using the formula cos(q + k) = cosgcos/c — sinqsinfc in Eq. (|15[) . we immediately see that the solution must be of 
the form f(k) = sin k/R E (k). A direct substitution then yields 



-2U r dq sm 2 q -2U f n dq sin 2 q 



E 2 J_ n 2nR E {q) 2 E 2 J_ n 2ir 1 - U/E - 2*7(1 + 2^ - t l cosq)/E 2 ' 



(16) 



where in the second equality we have made use of the asymptotic expansion of I E {q) from Eq.®. Next, we write 
E = U + a in Eq. ([TB| and take the limit U — > —oo assuming the energy shift a <C \U\. This yields 

/ n dq sin 2 q , . 

—L - (17) 

. 7r 2 7 r-a + 2(l + 2^)-2^cosg' v ' 
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FIG. 3: (Color online) Ratio between the trimer and pair energies as a function of t, calculated for decreasing values of the 
interaction U. Notice that for \U\ < 0.4 the curves are on top of each other signaling convergence to the continuum limit. We 
then recover the results for the continuum model obtained in Ref.— (see their Fig.l, left panel) with the (Bloch) masses given 
by m a = l/2t a . 



which can be solved analytically to obtain the shift a. By using the formula sin 2 q/(a — cosq)dq = 27r(a— V a 2 — 1) 
valid for a > 1, we find a — 4*4, — t 2 + 1. Making use of the strong coupling expansion E% T (U — > — oo) ~ — U — 2 — 2*j_, 
we finally obtain-^ 

£t(C/ = -oo) = (^-l) 2 , (18) 

showing explicitly that the binding energy of the trimer remains finite even in the strongly interacting regime. 

It is also interesting to consider the opposite limit of weak interaction, namely \U\ -C £4., if- In this case only 
the states at the bottom of the band are important and we can approximate the tight-binding dispersions with the 
quadratic expansions e CT (fc) ~ fc 2 /2?7v, where m a — l/2t <7 are the related Bloch masses. This corresponds to the 
continuum model studied in Ref^. The convergence to the continuum result is studied in Fig|3] where we plot the 
ratio \E\/Ep aiT between the trimer and pair energies as a function of *7 ^ 2 for decreasing values of the interaction 
strength \U\. We see that for the range of mass asymmetry considered here all curves with \U\ < 0.4 are on top of 
each other and coincide with the continuum prediction of Ref.—. Clearly deviations from the continuum limit are 
stronger for large mass asymmetry, where the condition — U <C t± becomes more stringent. 



B. Effective mass 



It is also interesting to discuss the effective mass M tI of the trimer. The latter is related to the energy dispersion 
E(P) of the trimer by 1/M tr = d 2 E/dP 2 evaluated at P — 0. We replace the derivative by a finite difference that 
we evaluate numerically. The result for the inverse effective mass is plotted in the inset of Fig. [4] as a function of 
the hopping rate tj. and for different values of the attraction strength. We see that the trimer becomes heavier as 
decreases or \U\ increases. 

At the symmetric point t± = 1, where the trimer disappears, the effective mass must coincide with the sum of 
the masses of its constituents. The effective mass of the pair can be calculated by the same formula, starting from 
the relation 1 + UIe{0) = 0, where the function Jg in Eq.© is evaluated at finite total quasi-momentum P ^ 0. 

This yields the energy dispersion E(P) = 2 + 2t^ — yV 2 + 4 + 4< 2 + 8t± cos P for the pair, from which we obtain 
= 4*4./ 4/4 + 8*4. + 4* 2 + U 2 . Since the effective mass of the heavy fermion is simply given by 1/(2*4.), we find 
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FIG. 4: Inverse effective mass of the trimer as a a function of the tunneling rate t± for different values of the interaction 
U — — 2(upper solid curve), — 4, — 8. The corresponding asymptotic curves [see Eq. (|19[l ] in the symmetric limit ij, = 1 are 
shown with the dashed lines. 



that the total mass M tot of the constituents is given by 

M tot = {^A^+Vf + U 2 + 2)/U u (19) 
which is shown in Fig |4] with dashed lines. 

IV. EXCITED BOUND STATES 

So far we have discussed the ground state solution corresponding to the bound state with the lowest energy (or, 
equivalently, with the largest binding energy). As mentioned in Section [lTIl Eq.(|9"j) admits other solutions corresponding 
to excited bound states, namely states with energy E — E n satisfying E\ < E„ < —E^ a ", where E — E\ is the ground 
state energy. 

In Figf5]we plot binding energy versus tunneling rate of the ground state (top curve) and the first four excited 
bound states calculated for U — —2. We see that these excited states are only stable for sufficiently large mass 
asymmetries, corresponding to tj, <C 1. In the limit = 0, Eq.® admits an infinite number of solutions of the form 
fn(k) — sin(rtfc). The corresponding energy levels E = E n are obtained from the condition 1 = y n (E), where the 
function y n (E) is defined in Eq.(JT5J). To understand this fact we notice that Rs(k) — Re is constant for t± = 0, so 
from Eq. (fTU| we find that 

1pn(ni,Tl2,ni) = — " [S n2 -m,-n - bn 2 - ni ,n] , (20) 

LiRe 

where A n is a normalization factor. Equation (I20p shows that the index n corresponds to the physical distance between 
the two (infinitely) heavy fermions. One is therefore left with the problem of a single light fermion in the presence 
of two static contact potentials separated by a distance d — n. An explicit solution shows that the ground state is 
bound for any distance d, and the corresponding energy is given by E — E n . It should be noticed that the case n = 0, 
corresponding to two heavy fermions at the same site, is forbidden in our three-body problem by the Pauli exclusion 
principle [the wave- function (|20|) vanishes]. As a consequence the ground state solution corresponds to n = 1, where 
the heavy particles are nearest neighbor. 

In the presence of a finite tunneling rates, tj, 0, the heavy fermions delocalize more and more affecting significantly 
the stability of these excited bound states, as shown in FigfSl 

Finally, it is interesting to discuss the behavior of the binding energy of the excited states as a function of the 
interaction strength. This is shown in FigJB]for a fixed value t± = 0.01 of the tunneling rates. We see that differently 
from the ground state solution, the binding energy for n > 1 exhibits a non-monotonic behavior as a function of U, 
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FIG. 5: (Color online) Binding energies of the ground state (upper curve) and the first four excited three-body bound states 
as a function of the hopping ratio and for U = —2. The excited bound states are only stable for large mass asymmetries 




\U\ 

FIG. 6: (Color online) Binding energies of the ground state (upper curve) and the first two excited three-body bound states 
as a function of the interaction \U\ for fixed t± = 0.01. Notice the non-monotonic dependence of the binding energy for the 
excited states n > 1. 



with a maximum around U ~ — 2. This comes from the fact that for \U\ 3> 1 the wavefunction for the light fermion 
is given by the superposition of two orbitals that are peaked at the positions of the heavy particles. When the latter 
are not nearest neighbor, corresponding to n > 1, the two orbitals have vanishing overlap as \U\ becomes large. As a 
consequence, the light fermion cannot easily delocalize between the two sites implying that the binding mechanism is 
less robust. 
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V. CONCLUSIONS 



In conclusion, we have presented a detailed discussion of three-body bound states of interacting fermions in one- 
dimensional optical lattice. These states only occur when the tunneling rates of the two spin components are different. 
As a consequence the asymmetric model (JTJ) exhibits a new strongly correlated phase at low but finite density, 
corresponding to a Luttinger liquid of trimers, which is completely absent in the usual attractive Hubbard model. This 
happens when the densities of the two components are commensurate allowing binding of all atoms into stable trimers. 
Then the interacting gas of trimers form a one-component Luttinger liquid. Thanks to the Pauli exclusion principle, 
these phases are expected to be particularly stable against three-body recombination and therefore experimentally 
accessible with ultra-cold Fermi gases. 

Our results provide one more example of the importance of few-body physics to identify the fundamental degrees 
of freedom of a full many-body system. 
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